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I.  INTRODUCTION 


In  the  development  of  a  predictive  tool  for  the  fluid  flow  field  due  to  the  interaction  of 
the  rocket  exhaust  plume  with  its  environment,  the  mixing  region  analysis  is  critical.  The 
manner  in  which  these  streams  interact  and  the  accurate  prediction  of  this  interaction  is 
paramount  to  several  missile  systems  applications.  Missile  signature  and  vehicle  design  are  two 
of  the  most  important  of  these  applications.  In  order  to  properly  simulate  the  mixing  region, 
the  numerical  calculational  procedure  must  be  accurate  and  the  turbulence  model  must  be 
physically  correct.  Having  a  physically  correct  turbulence  model  is  certainly  the  most 
demanding  of  these  two  requirements. 

Turbulence  has  been  investigated  for  many  years  now  and  is  usually  characterized  by  its 
randomness  and  disorderliness.  Despite  the  randomness  and  disorderliness  however, 
statistically  distinct  average  values  are  obtainable  for  the  velocity,  temperature,  and  density  for 
example.  The  randomness  and  disorderliness  is  characterized  by  scale  siz.e.  Not  only  is  the  fine 
scale  characterized  by  vortex  interactions  but  likewise  for  the  large  scale.  This  large  scale 
motion  has  been  studied  intensively  over  the  last  several  years  by  scientists  at  the  California 
Institute  of  Technology  and  offers  a  better  understanding  of  the  physical  phenomenology  of 
turbulence. 

II.  MATHEMATICAL  FLOW  MODEL 

The  mathematical  flow  model  utilized  in  this  investigation  consists  of  the  axisymmetric 
jet  mixing  equations  for  a  reacting  gas  mixture.  This  set  of  coupled  partial  differential 
equations  is  solved  utilizing  a  mixed  implicit/explicit  finite  difference  procedure.  The 
governing  equations  are  parabolic  and  are  solved  in  streamline  coordinates  using  a  marching 
scheme.  This  technology  was  essentially  developed  by  the  Joint  Army,  Navy,  NASA  and  Air 
Force  (JANNAF)  Plume  Technology  Working  Group  in  1972  with  the  development  of  the 
Low-Altitude  Plume  Program  (LAPP).  Technology  developments  since  that  time  have 
occurred  and  are  being  incorporated  in  the  JANNAF  Standardized  Plume  Flowfield  (SPF) 
Program.  The  improvements  to  the  mixing  portion  of  this  program  include  the  employment  of 
a  discretized  shear  layer  which  grows  due  to  the  mixing  of  thejet  and  the  external  streams.  This 
allows  a  more  optimum  placement  of  the  grid  points  in  the  flowfield.  Hence,  this  procedure  of 
retaining  the  ( x,  ib)  computational  grid  and  discretizing  the  shear  layer  has  led  to  a  much  more 
efficient  handling  of  the  numerical  procedures  used  to  solve  the  problem.  Other  improvements 
such  as  the  formulation  of  the  energy  equation  in  terms  of  total  enthalpy  rather  than  temp¬ 
erature  leads  to  more  accurate  solutions  in  higher  energy  rocket  propellants.  In  addition  a 
mass  flow  check  has  been  added  to  insure  that  mass  flow  is  truly  being  conserved.  This 
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model  lias  been  *'<  ::;undied  by  Aeronautical  Research  Associates  of  Princeton  (ARAP)and 
. ^iisiitiites  a  vital  working  portion  of  the  SP1’  program  being  developed  by  them  Cor  the 
JANNA1  Plume  Technology  Subcommittee.  This  eode  has  been  called  BOAT  in  previous 
references  in  the  literature  |9|.  Detailed  derivations  of  the  governing  fluid  dynamic  equa¬ 
tions  utilized  in  this  investigation  can  be  found  in  the  literature  [9  -  1  2 1  and  will  only 
briefly  be  presented  here.  The  governing  equations  arc 
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These  equations  are  then  transformed  from  the  (x.r)  to  the  (x,  ip)  coordinate  system  with  the 
transformation 
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Utilizing  the  transformation  in  (6).  the  governing  equations  become 
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I  vu>  turbulence  kinetic  energy  (TKF)  models  were  used  to  determine  the  turbulent 
viscosity  that  appears  in  the  governing  equations.  The  first  was  the  Kt2  model  developed  by 
Spalding  and  co-workers,  and  the  second  was  the  kai'  model  developed  by  Saffman  and  co¬ 
workers.  A  detailed  derivation  of  the  latter  model  development  for  this  investigation  is 
given  in  Section  IV  of  this  report.  Details  of  the  Kc2  model  can  be  found  in  other 
references  |  I  2 1 . 
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The  governing  turbulence  equations  for  the  ke2  model  are  given  by 
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Note  that  this  model  utilizes  five  empirical  constants  —  CV.  C«i,  CAj.  ak,  at. 

For  axisymmetric  flows  — 

(13) 
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When  these  corrections  are  made  to  the  constants,  the  model  is  known  in  the  literature 
as  the  kt  I  turbulence  model. 

It  was  determined  by  exercising  the  model  that  the  ke  1  model  could  not  accurately 
predict  weak  shear  flow,  i.e..  flow  in  which  the  two  streams  interacted  at  nearly  the  same 
velocities.  Therefore  a  correction  was  made  for  weak  shear  flows  by  altering  the  constant  C/i  as 
follows 
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NOTE:  f,  i  :  These  symbols  are  the  same 

throughout  report . 
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1  he  resulting  changes  in  the  ke  I  model  were  then  known  as  the  ke2  turbulence  model 
and  is  utilized  as  such  in  this  investigation. 


I  n  addition,  the  ke2  model  does  not  contain  any  terms  to  handle  compressibility  effects. 
Hence,  a  compressibility  correction  was  introduced  into  the  model  when  large  velocity 
differences  between  the  mixing  streams  became  important.  The  compressibility  correction  that 
was  used  resulted  from  an  empirical  formulation  due  to  Dash  1 1 1  ].  The  compressibility 
correction  factor  k  is  multiplied  by  the  constant  CM  and  ( 12)  becomes 
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where  k  is  a  function  of  the  maximum  turbulence  Mach  number 
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The  functional  form  of  k  is  shown  in  Figure  I. 


Figure  1.  Compressibility  correction  factor  for  kt2  turbulence  model. 
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The  k (o'  turbulence  model  was  also  utilized  in  this  investigation.  The  governing 
equations  for  this  model  are  developed  in  detail  in  Section  IV,  and  are  as  follows 
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This  model  also  uses  five  constants,  Cki,  Co.;,  C^i,  Gm.  Co-4,  and  when  compressibility 
effects  are  important,  an  additional  term  is  included  in  the  turbulence  kinetic  energy  equation 
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These  equations  were  transformed  to  the  (x,  \l>)  coordinate  system  utilizing  (6)  and  the 
following  equations  result  from  the  transformation  of  ( 10)  and  (II). 

For  the  ke2  model 

Turbulence  Kinetic  Energy  (TKE) 


3k 

3x 


r) 


(22) 
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(25) 


Similarly  lor  the  k«/  modeh  I XI and  (14)  we  transformed  lo  the(x.  iM coordinate  system 
with  the  result  that 


Turbulence  Kinetic  I  nergy  <  IK  I  ) 


Yk  =  1 
3x  Tty 


5  I  A  ak  I  + 

^  j  H  I 


U) 
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I  v  3V  | 
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e  p  1  e 
e 

Turbulence  I’seudo  Vorticity 
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where  A  is  delmed  in  (24| 


(26) 


(27) 


When  compressibility  ellects  are  not  important,  C is  taken  equal  lo  zero  Hence,  the 
compressibility  term  is  built  into  the  I  Kl  equation. 

Ill  STAR T1  INF  CONDITIONS 

I  hree  methods  were  utilized  to  define  the  startline  conditions  w  hich  are  used  to  initiate 
the  tinitc  hllerence  calculational  procedure  I  he  three  methods  are 
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Boundary  l  ayer  Initial  Profile 


•  Shear  1  aver  Initial  Profile 


Generalized  Specific  Profile 


I  he  first  two  met  hods  employ  calculational  procedures  to  generate  the  initial  profile.  The  third 
simply  specifies  the  value  of  all  the  variablesat  some  initial  downstream  location.  Each  of  these 
methods  was  utili/ed  to  calculate  the  shear  layer  flow  field  tor  the  two  dimensional  He  N:  case 
run  experimentally  bv  Brow  n  and  Roshko  |  1 1  for  a  density  ratio  of  1  -  and  a  velocity  ratio  ol  7. 

I  his  corresponds  to  the  experimental  data  given  in  figure  13a  of  reference  |  1 1. 


A  Boundary  l  ayer  Initial  Profile  Description 

I  he  displacement  effect  of  the  let  and  external  boundary  layers  is  calculated  by 
utilizing  a  velocity  profile  that  is  derived  from  the  combination  of  the  "l  aw  of  the  Wall"  and 
the  I  aw  of  the  Wake  "  hgure  2  shows  the  applicable  geometric  configuration. 


r  r.t- 


EXTERNAL  STREAM 


SPLITTER 

PLATE 


JET  STREAM 


*0  0 


Figure  2.  Boundary  layer  initial  profile. 


Note  that  the  shear  layer  starts  growing  at  x  =  Olrom  rc  =  r.  +6,  on  the  external  stream  side  and 
r  =  -6,  on  the  jet  stream  side.  The  profiles  utilized  to  calculate  these  initial  radii  are  computed 

from  the  velocity  piofiles 


=  -  2.5  In  £,  -  1.38  [2  -  w ( ' )  ] 

=  -  2.5  In  C  -  1. 38  [2  -  w(T)  ] 


(28) 


where  wi£i  is  C  ole's  universal  wave  function 


w(f.)  =  1  +  sin  n 


and  t  is  the  non-dimensional  radius 


C  =  |r  -  rc| 


e 

In  addition,  the  frictional  velocities 


(29) 


(30) 


Now  if  (he  displacement  thickness  is  known,  then  the  definition  of  this  quantity 


(32) 


15 


can  be  inverted  to  find  6,  provided  Ur  is  also  known.  Similarly  Sc  can  be  found  and  then  r,  and  rc 
can  be  evaluated.  Knowing  these  values,  the  initial  distribution  is  determinec 

B.  Shear  Layer  Initial  Profile  Description 

For  this  method,  a  fully  developed  shear  layer  profile  is  assumed  to  exist  at  the  initial 
streamwise  location.  The  initial  shear  layer  width  is  calculated  from  the  incompressible 
relation  |  1 2 1 


r 

e 


r  . 
J 


0.27 


u  .  + 
D 


(33) 


Figure  J  illustrates  this  notation 


Upon  establishment  of  the  shear  layer  upper  and  lower  boundaries  utilizing  (33).  the  properties 
are  distributed  across  the  shear  layer  according  to  the  simple  cubic  relations 


^j(K) 


Figure  3.  Shear  layer  initial  profile. 
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n  = 
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r 
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r  . 


r  . 
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I  he s he. 1 1  Li  v e i  grid  points  are  spaced  evenly  across  it  in  increments  of  At//  — 


h 


C.  tienerali/ed  Spectltc  Profile  Description 


I  his  method  ol  inputting  startline  conditions  does  not  rely  on  any  calculations  but 
meteh  uses  the  specified  profile  !  igurv  4  illustrates  this  case  where  u(  r),1\r),  and  x(r)are  input 
directly  at  the  startline  location  \  -  x„.  I  his  method  can  only  be  used  in  rare  cases  w'hen  a 
significant  amount  of  experimental  data  is  available. 


In  addition  to  these  fluid  dynamic  initial  conditions,  an  initial  turbulence  level  must  be 
supplied  ai  the  initial  axial  station  \...  In  the  absence  of  known  profiles  for  k  and  f  or  u' 
the  I’randil  mixing  length  model  is  used  to  define  the  turbulent  shear  stress  in  terms  of  the 
local  velocity  gradient  utilizing  the  following  relation 


4 


I 


p  U  V 


3u 

3r 


(35) 


'  o2  /3u  \ 

-  (37) 


and  the  relation  between  the  shear  stress  and  the  turbulent  energy  [7] 


1  '  ' 

k  =  u  v 


0.3 


(36) 


D.  Boundary  Layer  Initial  Profile  Results 

Lor  the  He  N;  shear  layer  of  Brown  &  Roshko,  the  initial  conditions  should  not  affect 
the  resulting  similar  profile  far  downstream  of  the  initial  station.  This  was  examined  by 
looking  at  the  different  methods  of  specifying  the  initial  profile. 

The  use  of  the  boundary  layer  initial  profile  was  examined  for  the  He/N:  shear  layer  of 
Brown  and  Roshko  for  which  the  experimental  velocity  and  density  profiles  are  shown  in 
Figures  lb  anil  I  7. 

6*  was  determined  from  6  calculated  by  Brown  and  Roshko  and  the  relation 
between  these  quantities  for  a  flat  plate,  i.e. 

6*  _  1.7208  Re^  =  2  592 
6  0.664  Re15 


j 


Since  the  momentum  thickness  was  found  to  be  0.001  inches  and  the  splitter  plate  thickness 
was  0.002  inches,  the  displacement  thickness  utilized  for  this  option  was 

6*  =  (0.002)  (2.592)  =  4.318  x  10~4ft. 

- u 

The  Re  difference  between  the  He  and  N*  was  not  accounted  for  and  the  boundary  layer 
displacement  thickness  was  taken  as  the  same  for  both  the  jet  and  the  external  stream,  i.e. 

<5*  =  8*  =  4. 32  x  10-4ft. 

3  e 

The  resulting  mixing  profiles  of  the  density  and  the  velocity  for  the  two  dimensional 
shear  laser  arc  shown  in  Figure  6.  These  results  are  given  at  a  distance  downstream  of  x  =  1.69 
inches. 
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Figure  5.  Results  for  boundary  layer  initial  profile. 


I  Shear  l  ayer  Initial  Profiles  Results 

l  h is  method  was  used  to  ealeulate  the  initial  profile  from  which  the  finite  difference 
solution  was  started  and  the  results  from  this  method  are  shown  in  Figure  7.  The  initial  shear 
laser  based  on  the  cubic  profile  was  calculated  at  x  —  .01  feet  (8.33  x  10  J  inches)  and 
calculations  were  started  from  that  point.  The  kt2  turbulence  model  was  used  in  this 
calculation  just  as  it  was  on  the  boundary  layer  initialization  results  given  in  Section  Dabove. 
Note  that  these  results  are  given  for  a  position  slightly  further  downstream  x  =  1.81  inches. 
Ibis  makes  very  little  difference  however,  since  the  profiles  have  already  achieved  self- 
smnlai  it\ 

K.  Generalized  Specific  Profile  Results 

This  method  was  also  used  to  specify  the  initial  profile  to  determine  what  differences  in 
results  occurred  because  of  differences  in  the  initial  startline.  For  this  option,  the  mean  profile 
for  the  density  and  velocity  were  taken  at  x  -  0  directly  from  the  boundary  layer  initialization 
scheme  described  above.  However,  the  turbulence  kinetic  energy  profile  at  x  =  0  was  taken 
from  experimental  data  rather  than  being  calculated  using  the  mixing  length  turbulence 
option  described  above. 

I  he  turbulence  kinetic  energy  profile  was  taken  from  experimental  data  for  flow  overa 
flat  plate.  This  data  was  taken  from  Figure  8  of  Klcbanoff  1 3 1  and  has  been  tabularized  in  Table 
1  as  a  function  of  v.'fi  where  8  is  the  shear  laver  width. 


TABLE  1.  EXPERIMENTAL  TURBULENT  KINETIC  ENERGY  PROFILE 


jTiiYLusdli 

U  2 


1.2818x10-2  6.409x10-3 

9.345x10-3  4.673x10-3 

8.456x10-3  4.228x10-3 

7.348x10-3  3.674x10-3 

6.548x10-3  3.274x10-3 

5.371x10-3  2.686x10-3 

3.830x10-3  1.915x10-3 

2.481x10-3  1.241x10-3 

1.206x10-3  6.030x10-4 

4.320x10-4  2.160x10-4 

1.470x10-4  7.350x10-5 


Next,  the  velocity  profile  that  was  generated  at  x  =  0  using  the  boundary  layer 
initialization  scheme  above  is  interpolated  to  obtain  velocities  at  each  y/<5  location.  This  was 
done  for  both  the  jet  boundary  and  the  external  boundary.  These  results  are  shown  in  Tahir  2. 
Also  note  that  the  initial  shear  layer  thickness  6  is  taken  from  the  boundary  layer  initialization 
scheme.  The  combination  of  the  data  for  both  the  external  and  jet  streams  provide  the  specific 
shear  layer  profile  data  used  to  calculate  the  mixing  layer.  This  is  given  in  Table  3. 

This  is  compared  with  the  initial  turbulence  data  which  was  calculated  utilizing  the 
mixing  length  model  as  described  above.  This  initial  profile  expanded  to  50  points  across  the 
shear  layer  is  shown  for  the  boundary  layer  initialization  method  in  Table  4  and  for  the  specific 
profile  method  in  Table  5. 

Comparing  Table  4  with  Table  5  ..  ,ows  large  differences  in  the  turbulence  kinetic 
energies  for  the  calculated  and  input  initialization  schemes.  Note  that  for  the  calculated 
scheme  of  Table  4,  0  sS  k  ^  792.6  while  the  data  based  on  experiment  is  in  the  range  0.079^ 
5.657.  Hence,  the  turbulent  intensity  is  down  two  orders  of  magnitude.  Similarly,  there  are 
large  changes  noted  in  length  scale  parameter  e.  However,  when  these  shear  layers  have  been 
calculated  out  to  a  distance  of  1 .69  inches.  Figure  (7) shows  that  the  initial  profile  differences 
have  washed  out  and  the  density  and  velocity  profiles  are  virtually  the  same.  Figure  (H)  was 
obtained  from  Figures  (5-7). 

Therefore  it  has  been  shown  that  this  calculational  scheme  is  not  sensitive  to  initial 
conditions  when  the  behavior  of  the  shear  layer  isexamined  far  enough  dow  nstream  where  the 
How  becomes  self  similar.  Hence,  any  of  the  initial  profile  methods  may  be  used  with 
confidence. 


IV.  k to'  TURBULENCE  MODF.I  FORMULATION 

In  order  to  make  meaningful  comparisons  of  various  turbulence  models  lor  use  in 
rocket  exhaust  plumes,  an  investigation  of  several  models  was  made.  The  following  turbulence 
models  were  investigated: 

•  I’randtl  mixing  length 

•  Donaldson-(iray  eddy  viscosity 


23 


r 


TABLE  2.  CALCULATION  OF  TURBULENCE  KINETIC  ENERGY 
ACROSS  JET  AND  EXTERNAL  BOUNDARIES  USING 
EXPERIMENTAL  DATA  FROM  TABLE  1 


JET  BOUNDARY 


y/fi 

k/U^ 

Use 

yt 

k 

ut 

1- 

0.0 

6.409x10-3 

10.00000 

6.8992 

0.0 

300 

0.1 

4.673x10-3 

9.999847 

5.0305 

25.84 

B9 

4.228x10-3 

32.81 

9.999694 

4.5514 

27.462 

300 

E3 

3.674x10-3 

32.81 

9.999541 

3.9550 

28.541 

300 

0.4 

3.274x10-3 

32.81 

9.999388 

3.5244 

29.449 

300 

0.5 

2.686x10-3 

32.81 

9.999235 

2.8914 

30.255 

300 

0.6 

1.915x10-3 

32.81 

9.999082 

2.0615 

30.980 

300 

0.7 

1.241x10-3 

32.81 

9.998929 

1.3359 

31.610 

300 

0.8 

6.030x10-4 

32.81 

9.998776 

0.6491 

32.136 

300 

0.9 

2.160x10-4 

32.81 

9.993623 

0.2325 

32.536 

300 

1.0 

7.350x10-5 

32.81 

9.99847 

0.0791 

32.810 

300 

EXTERNAL  BOUNDARY 


y/<5 

k/U^ 

Uoc 

y 

k 

ut 

T 

0.0 

6.409x10-3 

4.69 

10.00000 

0.14097 

0.0 

300 

0.1 

4.673x10-3 

4.69 

10.00015 

0.10279 

3.6915 

EM 

4.228x10-3 

4.69 

10.00030 

0.09300 

3.9241 

300 

KB 

3.674x10-3 

4.69 

10.00045 

0.08081 

4.0800 

300 

0.4 

3.274x10-3 

4.69 

10.00060 

0.07202 

4.2094 

300 

0.5 

2.686x10-3 

4.69 

10.00075 

0.05908 

4.3247 

300 

0.6 

1.915x10-3 

4.69 

10.00090 

0.04212 

4.4283 

300 

0.7 

1.241x10-3 

4.69 

10.00105 

0.02730 

4.5186 

300 

0.8 

6.030x10-4 

4.69 

10.00120 

0.01326 

4.5933 

300 

0.9 

2.160x10-4 

4.69 

10.00135 

0.00475 

4.6508 

300 

1.0 

7.350x10-5 

4.69 

10.00150 

0.00162 

4.6900 

300 

tFrom  the  He/^  B.L.  initialization  ouput- 

6jet  =  10.0000  -  9.9985  =  1.5x10'3ft  (0.018  inches) 
5ext  =  10.0020  -  10.0000  =  2.0x10'3ft  (0.024  inches) 

UNITS:  y  -  in.;  U  -  ft/sec;  k  -  ft2/sec2;  T  -  °K. 
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TABLE  3. 


SHEAR  LAYER  INPUT  PROFILE  DATA  FOR  BROWN 
AND  ROSHKO  He/N2  EXPERIMENTAL  RUN  (FIGURE  16) 


PT. 

y 

; 

T 

k 

"He 

"N2 

1 

9.998470 

32.810 

300.00 

0.0791 

1.0 

0.0 

2 

9.998623 

32.536 

300.00 

0.2325 

1.0 

0.0 

3 

9.993776 

32.136 

300.00 

0.6491 

1.0 

0.0 

4 

9.993929 

31.610 

300.00 

1.3359 

1.0 

nnm 

5 

9.999082 

30.980 

300.00 

2.0615 

1.0 

6 

9.999235 

30.255 

300.00 

2.8914 

1.0 

Jtffgps 

7 

9.999388 

29.449 

300.00 

3.5244 

1.0 

8 

9.999541 

28.541 

300.00 

3.9550 

1.0 

KIM 

9 

9.999694 

27.462 

300.00 

4.5514 

1.0 

pffflj 

10 

9.999847 

25.840 

300.00 

5.0305 

1.0 

HlSti 

11 

10.000000 

0.000 

300.00 

68992 

0.5 

12 

10.000150 

4.6900 

300.00 

0.10279 

0.0 

1.0 

13 

10.000300 

4.6508 

300.00 

0.09300 

0.0 

1.0 

14 

10.000450 

4.5933 

300.00 

0.08081 

0.0 

1.0 

15 

10.000600 

4  5186 

300.00 

0.07020 

0.0 

1.0 

16 

10.000750 

4.4283 

300.00 

0.05908 

0.0 

1.0 

17 

10.000900 

4.3247 

300.00 

0.04212 

0.0 

1.0 

18 

10.001050 

4.2094 

300.00 

0.02730 

0.0 

1.0 

19 

10.001200 

4.0800 

300.00 

0.01326 

0.0 

1.0 

20 

10.001350 

3.9241 

300.00 

0.00475 

0.0 

1.0 

21 

10.001500 

3.6915 

300.00 

0.00162 

0.0 

1.0 

UNITS:  y  —  in.;  U  -  ft/sec;  k  -  ft^/sec^;  T  -  °K. 
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TABLE  4.  INITIAL  PROFILE  FOR  BROWN  AND  ROSHKO  He/No 
SHEAR  LAYER  (FIGURE  16 )  —  BOUNDARY  LAYER 
INITIALIZATION 
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•  I  aundei -Spalding  kr2  turbulence  kinetic  cncigv  model 


•  Sallman  k<«  turbulence  kinetic  cnci p\  model 

I  lie  Sallman  model  was  chosen  because  this  model  has  been  formulated  to  account  lor 
a  speed ic  physical  phenomena  which  is  not  explicitly  modeled  in  the  other  three  1  his  model 
was  formulated,  primarily  on  empirical  arguments,  to  account  lot  shear  Mows  which  do  not 
have  a  constant  density  I  hese  density  differences  arise  either  due  to  differences  in  molecular 
weights  of  the  shear  laser  tlutds  (  heterogeneous  lluidsi  or  due  to  compressibility  effects  (Mach 
Number)  in  homogeneous  fluids  Both  ol  these  effects  are  modeled  in  the  turbulence  equations 
formulation 

compressibility  etfects  are  accounted  lor  in  the  kr2  turbulence  model,  but  in  an  ad  hoc 
manner  I  he  turbulent  viscosity  contains  an  empirical  correction  term  k(M.,„,,i  which  was 
formulated  using  only  a  very  limited  set  of  experimental  data  1  here  are  no  additional  terms  m 
the  model  that  account  lor  this  effect. 

In  order  to  use  the  Sallman  model  to  make  comparisons  with  data  it  was  necessary  to 
extend  this  model  to  a  cylindrical  geometry  and  to  formulate  it  using  a  t mite  dillerencing 
scheme  common  to  the  other  three  models  investigated  Only  in  this  wav  can  differences  in 
results  be  solely  due  to  the  turbulence  modeling  and  not  due  to  the  numerics.  In  addition,  the 
equations  are  formulated  in  a  stream  function  coordinate  system 

C  onsider  the  Sallman  formulation 
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w  here 


A  =  0.09  a"  =  B'=  5/3  a'  =  B'a"  -  2k2 

I  liese  equations  can  he  generalized  for  either  2-D  or  axisymmetric  flow  to  be 
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In  order  to  apply  these  equations  to  an  axisymmetric  roc  I ’t  exhaust  plume,  choose  j  -  1: 
y  -  r.  L  tili/ing  different  notations  tor  the  turbulence  kinetic  energy  ,k  ~  11 '  *  .1*  '•  I*lcsc 
equations  become 
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where  c  k  &  u  = 
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so  rewriting  these  equations. 
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[  xpanding  ( I )  for  a  steady  flow  gives 


This  equation  can  be  transformed  from  (x.r)  —  utilizing  a  transformation  of  dependent 
\ariables. 
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Expanding  (38)  for  steady  flow  gives  — 
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where: 
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and  the  other  constants  are  as  they  appear  on  page  34.  To  account  for  the  effects  of 
compressibility.  Saffman  introduces  an  additional  term  into  the  k-equation.  This  is  the  last 
term  in  his  equation  given  below 


where  the  subscript  “1"  in  the  last  term  refers  to  properties  of  the  external  stream. 
Rewriting  this  equation 
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Substituting  this  into  the  previous  equation  gives 
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Now  if  we  utilize  the  notation  that  the  turbulent  kinetic  energy  is  given  by  k,  Equation  (48) 
becomes 
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where  as  before  k  =  e 

writing  this  equation  in  terms  of  the  substantial  derivative 
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Now  for  a  perfect  gas,  the  sonic  velocity  of  the  external  stream,  ai  is  given  by 


Also  for  a  perfect  gas 
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Where  the  subscript  (  refers  to  the  externa!  stream.  Utilizing  this  relation  (51)  in  (50)  above 
leads  to 
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Transforming  this  equation  to  cylindrical  coordinates 
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Now  if  we  define 
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and  drop  the  mean  value  notation  (53)  becomes 


(55) 


where 


Apk  _  Apk 


(56) 


Utilizing  (56)  the  second  term  on  the  right 


hand  side  (RHS)  of  (55)  becomes 
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Substituting  (57)  back  into  (56) 
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Equation  (58)  becomes 
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Where  (59)  differs  from  (47)  in  the  previous  formulation  only  by  the  last  term. 
Expanding  (59)  for  a  steady  flow  gives 
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Now  (60)  can  be  transformed  from  (x.r)  —  (x,</0  using  the  transformation  of  dependent 
variable 
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(62) 


(63) 
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Since  pu  always  >  0  for  the  shear  flows  considered  here 
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Where  it  will  be  noted  that  (65)  differs  from  (54)  only  in  the  addition  of  the  last  term 
From  Equations  (65)  and  (56)  — 
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where  the  subscripts  “2”  and  “ex”  are  synonomous.  Now  let  A  =  (pt+pur2)/^  and  utilize 
the  following  finite  difference  formulation 


42 


Hence  equations  ( h9  land  i  72  i  are  utilized  to  calculate  the  changes  in  k  and  /  along  the 
marching  direction  at  mesh  points  inside  the  calculational  lie  Id  Howcvei.  special  treatment  ol 
these  equations  is  neeessarv  along  the  axis  ol  the  Mow  as  shown  below 
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Now  along  the  axis  r  =  tg  =  0 
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Hence,  the  first  and  second  terms  on  the  RHS  ol  i73i  must  he  evaluated  along  the  axis  since 
these  terms  are  indelimte 
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Apple  ing  I  Hospital  s  rule  to  a  portion  ol  this  term  gives 
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I  Ik  lust  two  terms  on  i lie  R  US  ol  ( 75 1  arc  analogous  to  the  first  two  terms  on  the  R  HS  of  (74), 
ami  are  evaluated  analogously  lienee  looking  at  the  third  term 
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The  fourth  term  remains  as  is  and  the  /-equation  along  the  axis  becomes  - 
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Along  the  axis  (79) 


Hence,  evaluation  of  the  turbulence  model  equations  for  the  case  when  r  =  o  along  the  axis 
becomes  according  to  (77)  and  (7X) 


Now  when  compressibility  effects  are  important  the  last  term  in  the  k-equation  must  also  be 
evaluated  along  the  axis 
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Investigate  this  term  in  the  limit  as 
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Rew  rite  <  SO )  as 
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I  hen  using  I. ’Hospital's  rule  to  (XI) 
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Hence  in  the  limit  as  r  —  o  this  term  adds  nothing. 

Therefore  equations  (69)  and  (7 2)  lor  field  mesh  points  and  (84)  and  (85)  lor  points 
along  the  axis  determine  the  turbulence  kinetic  energy  and  the  pseudo  vorticity  in  the  shear 
laser. 


T  his  formulation  was  coded  and  added  to  the  BOAT  portion  of  the  SP1  code  now 
under  development.  These  coding  changes  were  input  via  an  update  to  the  main  code  and  are 
detailed  in  Appendix  B. 
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Along  the  axis 


V.  NON-RI  ACTING  SH MAR  l.AYF.R  COMPARISON 

In  order  to  evaluate  the  various  turbulence  mixing  models,  predictions  were  made 
corresponding  to  the  careful  experimental  measurements  made  by  Brown  and  Roshko 
|  I  |  over  an  extended  period  of  tune  The  experiments  were  made  in  a  laboratory  utilizing 
a  splittei  plate  separating  two  4  x  1-inch  2-0  nozzles  1  he  principal  aim  of  this  experimental 
work  was  to  investigate  the  effect  ol  density  differences  between  the  two  mixing  layers.  This 
was  accomplished  experimentally  by  using  various  mixtures  of  He  and  V.  Ar  was  also  used  in 

U  x 

some  rare  instances.  These  two  streams  were  turbulent  with  the  Re  number  — ^ —  in  the 

range  of  c  10'’.  The  experimental  device  was  run  at  low  speeds  in  the  range  of  «=  50  fps. 
freest  ream  velocity  and  density  ratios  on  the  order  of5®  7  were  run  experimentally.  (his  work 
showed  that  the  large  structure  existed  over  all  the  density  ranges  tested. 

Predictions  were  made  for  all  the  experimental  runs  made  by  Brown  and  Roshko  and 
are  presented  subsequently  Velocity  and  density  profiles  were  calculated  as  a  function  of  y  (x- 
x., I  where  x  is  the  virtual  origin  of  the  shear  layer  and  y  is  the  distance  above  or  below  the 
dividing  streamline.  The  dividing  streamline  was  located  utilizing  a  numerical  integration 
scheme  I  he  details  and  limitations  ol  this  calculation  is  given  in  Appendix  A. 

In  addition  to  the  velocity  and  density  profiles  that  were  compared  for  the  data  of 
Brown  anil  Roshko.  the  growth  of  the  shear  layei  asa  function  of  velocity  ratio  was  compared. 
As  a  basis  lor  this  comparison,  they  used  the  velocity  profile  maximum  slope  thickness  and  its 
derivative 
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Filiation  (87l  was  used  for  the  comparison. 

Results  comparing  the  shear  layer  flow  utilizing  the  two  turbulence  kinetic  energy 
models  with  experimental  data  are  shown  in  f  igure  V  I  his  is  a  comparison  of  the  constant 
density  air  air  case  w  it h  a  jet  velocity  of  32.8  feet  second  and  an  external  stream  \  .locity  of  4.7 
feet  second.  Ihe  flow  was  at  a  constant  pressure  of  102.9  psia. 

Note  that  the  agreement  between  theory  and  experiment  is  very  good  for  both  the  kt2 
and  the  kcc'  turbulence  models.  I  he  shear  width  agreement  is  excellent  and  velocity  profile 
slope  is  very'  good.  The  kt2  model  shows  a  little  belter  agreement  with  the  slope  while  the 
kcc'  model  shows  a  slightly  better  shear  width  agreement.  Overall  it  can  be  said  that  the 
agreement  with  both  is  very  good. 

I  he  next  comparison  made  was  another  constant  density  air  air  case  for  which  the 
pressure  was  maintained  constant  at  102.9  psia  For  this  case,  the  jet  velocity  was  held  at  32.8 
feet  second  while  the  external  stream  velocity  was  increased  to  12.5  feet  second.  This 
comparison  is  shown  in  Figure  10. 

Note  that  the  agreement  between  theory  and  experiment  is  again  pretty  good  for  both 
turbulence  models.  Ihe  slope  is  reasonable  for  both  and  the  width  of  the  shear  layer  is 
approximately  the  same  for  both.  The  ku/  model  falls  closer  to  the  actual  data  than  does  the 
kt2  model. 

Hence,  for  the  constant  density  turbulent  2-1)  shear  layer,  the  turbulence  prediction 
models  are  doing  reasonably  well  at  predicting  the  growth  and  velocity  profiles. 

In  order  to  evaluate  shear  layer  flows,  where  multiple  species  are  involved  such  as  for 
the  lie  \  experiments,  it  is  necessary  to  know  the  value  of  the  turbulent  I’randtl  number.  This 
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Velocity  profile  comparison  for  air/air  shear  layer 
Table  6,  case  number  I. 
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Velocity  profile  comparison  for  air/air  shear  layer  -  Table  6 
case  number  II. 


parameter  arises  in  the  governing  differential  equations  to  account  lor  differences  belweer  the 
velocity  profile  and  the  specie  profiles.  I  his  parameter  is  not  known  a  priori.  It  can  be 
determined  however  from  the  experimental  data. 


Brown  and  Roshio  1 1 1  found  that  for  all  cases  of  Hc/N,  shear  layers,  the  spreading 
angle  of  the  density  profile  was  greater  than  that  for  the  velocity  profile.  In  order  to  evalu¬ 
ate  this,  they  constructed  an  eddy-viscosity  model  and  deduced  that  the  turbulent  Prandtl 
number  should  be  between  0.2  and  0.3. 

l-ixurv  II  compares  the  predicted  pu  profile  for  the  conditions  of  Case  V  shown  in 
I'ablc  n  utilizing  the  present  kt2  turbulence  kinetic  energy  model.  Note  that  there  is  qualitative 
agreement  between  the  predictions  and  the  experiment  for  a  turbulent  Prandtl  number  of  0.3. 
Had  this  been  for  the  constant  density  case,  the  functional  relationship  would  have  been  much 
different,  approaching  a  constant  value  of  !.().  The  trend  toward  matching  the  experimental 
data  is  to  run  at  lower  Pr,  numbers.  However,  this  studs  did  not  investigate  the  quantitative 
differences  as  a  lower  Pr  is  utilized.  One  reason  for  this  was  the  marked  increase  in  computer 
runtime  that  would  have  been  necessary.  This  is.  however,  a  parameter  that  needs  to  be 
investigated  in  future  investigations  of  turbulence  modeling.  Therefore.  Pr,  =  0.3  was  chosen 
for  the  corresponding  calculation. 


TABLE  6.  INITIAL  CONDITIONS  FOR  SHEAR  LAYER 
COMPARISON  CASES 


CASE 

NO. 

UjAJe 

pj/pe 

JET  STREAM 

CONSTITUENT 

EXTERNAL 

STREAM 

CONSTITUENT 

1 

7.00 

1.00 

AIR 

AIR 

II 

2.62 

1.00 

AIR 

AIR 

III 

2.65 

7.00 

n2 

He 

IV 

2.65 

0.143 

He 

n2 

V 

7.00 

0.143 

He 

n2 

A  comparison  of  the  two  turbulent  kinetic  energy  models  with  experiment  is  shown  in 
terms  ot  the  velocity  profile  in  figure  12  tor  case  111  (  Table  6).  For  this  case,  the  jet  velocity  was 
32.8  feet  second  while  the  external  stream  velocity  was  12.4  feet  second.  The  jet-fluid  was  N: 
while  the  external  lluid  was  Fie  giving  a  density  ratio  of  7.  As  before,  the  pressure  was  held 
constant  at  1 02. *3  psia. 

Both  models  show  reasonable  agreement  with  the  data.  The  velocity  profile  slope  is 
more  nearly  constant  tor  the  kt2  as  opposed  to  the  ku/  model.  There  is  some  disagreement  with 
the  width  of  the  mixing  layer  between  the  models  where  the  width  is  too  great  on  the  high 
velocity  side  on  one  and  too  narrow  on  the  low  velocity  side  on  the  other  and  vice  versa. 

The  first  density  profile  comparison  is  shown  in  Figure  13.  (Case  111)  The  most  notable 
aspect  of  this  comparison  is  the  lack  of  agreement  between  theory  and  experiment.  In 
particular  the  slope  of  density  profile  on  the  N:  side  is  far  too  large:  the  kt2  model 
demonstrating  the  worst  agreement  between  the  two.  The  width  of  the  density  layer  is  very 
close  for  the  ku/  model  and  somewhat  worse  for  the  kr2  model.  The  absolute  agreement 
between  experiment  and  theory  is  poor  everywhere  across  the  mixing  layer  and  density  errors 
of  ~  100  percent  can  be  seen.  Flad  the  disagreement  occurred  only  in  the  edge  regions  of  the 
shear  layer,  concern  lor  this  would  have  been  lessened.  Unfortunately,  the  agreement  is 
uniformly  poor. 

The  next  comparison  of  the  turbulent  kinetic  energy  models  with  experiment  is  show  n 
in  terms  of  the  velocity  profile  in  Figure  14  for  case  IV  (  Table  6).  For  this  case,  the  jet  velocity 
was  32.8  feet  second  while  the  external  stream  velocity  was  12.4  feet  second.  The  jet  fluid  was 
Fie  w  hile  the  external  fluid  was  N;  giving  a  density  ratio  of  0. 143.  the  external  stream  being  the 
more  dense.  Again  the  pressure  was  held  constant  at  102.9  psia.  A  turbulent  Prandtl  number  of 
0.7  was  used  in  the  predictions. 

For  this  case,  the  velocity  profile  comparison  begins  to  look  poor,  especially  on  the  Fie 
side  of  the  shear  layer.  Some  of  this  poor  agreement  can  be  attributed  to  being  in  the  edge 
region  ot  the  shear  layer.  However,  it  is  clear  that  this  is  not  the  only  reason  for  the 
disagieement.  Notice  that  the  shear  layer  width  is  predicted  much  more  narrow  than  that 
observed  experimentally.  Further,  the  velocity  profile  slope  on  the  He  side  is  considerably  in 
error.  This  velocity  profile  comparison  is  the  worst  that  has  occurred  so  far. 

The  density  profile  comparison  for  C  ase  IV  is  show  n  in  Figure  15.  The  same  problems 
that  were  evident  in  the  velocity  profile  comparison  arc  magnified  for  the  density  profile.  The 
predicted  density  width  is  too  narrow  and  discrepancies  are  most  notable  in  the  slope  of  the 
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Velocity  profile  comparison  for  He/N2  shear  layer,  Table  6  - 
case  number  IV. 


density  profile  on  the  V  side.  The  differences  in  width  between  the  velocity  and  density  profiles 
was  accounted  for  in  the  model  by  running  the  model  at  a  turbulent  Prandtl  number,  Pr,=0.3 
as  described  earlier.  It  is  obvious  that  the  resulting  theoretical  dilference  in  width  of  the  two 
lasers  is  tar  less  than  the  experimental  width  dilference.  Thus,  it  is  obvious  that  there  are  some 
serious  problems  in  the  turbulence  modeling  for  flows  have  a  large  density  dilference. 

|  he  next  comparison  of  the  turbulent  kinetic  energy  models  with  experiment  are  the 
most  interesting  ol  the  shear  flows  compared  since  the  density  ratio  of  7  that  was  run  is  very 
nearly  that  seen  in  rocket  exhaust  plume  firings.  I  he  velocity  protile  is  shown  in  Figure  16  for 
Case  V  i  Fable  6).  for  this  case,  the  jet  velocity  was  32.8  feet  second  while  the  external  stream 
velocity  was  4.7  feet  second.  The  jet  fluid  was  He  while  the  external  lluid  was  N:  giving  a 
densitv  ratio  of  7,  the  jet  stream  being  the  less  dense  stream.  The  pressure  was  held  constant  at 
102.9  psia. 

1  he  velocity  profile  comparison  for  this  case  is  show  n  in  Figure  16.  As  was  the  situation 
in  Case  IV.  the  velocity  profile  comparison  looks  poor,  especially  on  the  He  side  of  the  mixing 
layer.  Note  the  substantial  difference  between  the  velocity  profile  slopes  predicted  by  both 
models  and  measured.  Again  the  shear  layer  width  predicted  is  more  narrow  than  that 
measured.  The  comparison  is  very  similar  to  that  of  Case  IV. 


The  density  profile  comparison  for  Case  V  is  shown  in  Figure  17.  Again  as  in  Case  IV. 
the  predicted  density  width  is  too  narrow  and  discrepancies  are  most  notable  in  the  slope  of  the 
density  profile  on  the  N;  side.  The  differences  in  width  between  the  v  elocity  and  density  profiles 
were  accounted  for  in  the  model  by  running  at  a  Pr,=0.3.  Again  the  predicted  width  difference 
between  the  predicted  density  and  velocity  profiles  are  far  less  than  between  the  measured 
profiles.  The  comparison  between  experiment  and  theory  is  the  worst  thus  far  seen  and  as 
mentioned  earlier,  this  is  the  case  of  most  interest  since  it  more  closely  matches  a  real  rocket 
plume  in  terms  ol  density  ratio. 


I  he  spreading  rate  for  all  the  previous  cases  was  calculated  and  compared  with  the 
experimental  data.  Figure  IS  shows  the  comparison  for  the  kt2  model  and  the  comparison  for 
the  k«/  model  is  shown  in  Figure  19.  There  is  one  difference  between  the  velocity  profile  slope 

(~\ 

of  the  experimental  data  and  the  theoretical  calculations  Brown  and  Roshko  used\<5y  /M  AX 


/mi\ 


/bu\ 


while  V>v/\V(i  was  utilized  for  the  theoretical  calculations.  In  calculating  the  \dy/AV(i  . 
onlv  the  center  50  percent  of  the  profile  was  utilized  in  determining  the  average  slope  to 
minimize  edge  gradient  ellecls. 
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It  should  be  noted  that  the  kt2  model  gave  the  best  agreement  between  theory  and 
experiment.  I  or  the  ease  with  equal  densities  lor  both  the  jet  and  external  streams,  the  data 
points  tall  exactly  on  the  theoretical  curve. 

The  agreement  between  theory  and  experiment  is  the  worst  for  the  case  where 
p,  n  -  1  ■  for  both  models. 

I  he  agreement  or  lack  thereof  between  experiment  and  theory  is  almost  totally 
governed  b>  jj~  .  Since  the  spreading  rate  is  so  sensitive  to  this  parameter,  these 
comparisons  are  much  less  meaningful  than  the  density  profile  for  example. 

VI.  NON-REACTING  JET  COMPARISON 

I  lie  previous  comparisons  of  the  turbulence  mixing  models  were  made  utilizing  the 
experimental  data  of  Brow  n  and  Rosliko  which  were  for  two  dimensional  shear  layers  at  low 
velocities  (.10  feet  second  or  less).  Since  the  applications  of  interest  for  this  work  are  all  at 
much  higher  velocities  and  since  the  geometry  is  axisymmetric.  it  was  felt  that  comparison  with 
some  of  the  NASA  Shear  Mow  Conference  Data  |7|  was  in  order.  Hence  comparisons  were 
made  for  axisymmetric  jet  data  in  order  to  compare  the  turbulence  models.  Two  sets  of 
experimental  data  were  chosen  from  the  NASA  Shear  Flow  Conference  for  comparison  with 
the  two  turbulence  kinetic  energy  models,  the  kt  2  and  k<y'.  Table  7  details  the  flow  conditions 
that  were  run  during  the  experiments. 

These  two  cases  cover  the  spectrum  of  expected  velocity  and  density  ratios  that  one 
might  expect  to  see  in  a  realistic  rocket  plume  case.  They  are.  however,  not  in  the  same 
experiment.  For  both  of  these  cases  only  experimental  velocity  profiles  were  measured.  This  is 


TABLE  7.  INITIAL  CONDITIONS  FOR  JET  MIXING  COMPARISON  CASES 


CASE 

NO. 

=: 

3  • 

p  .  /  p 

j  e 

JET  STREAM 
CONSTITUENT 

JET 

STREAM 

CONSTITUENT 

.  .  .  _ 

M  . 

J 

M 

e 

I 

no 

1  .97 

A  1  R 

A I  R 

2  .  2 

0 

I  I 

. 

2  .  72 

0.06 

H  2 

AIR 

0 . 89 

1  .  32 

6‘) 


unfortunate  since,  as  was  seen  in  the  preceding  comparisons  tor  the  shear  layer,  the  density 
profiles  are  a  much  more  stringent  test  for  the  accuracy  ol  the  theoretical  models, 
furthermore,  since  species  concentration  is  the  one  ol  the  quantities  used  directly  in  rocket 
plume  applications,  it  is  a  more  important  measure  of  the  accuracy  of  the  turbulence  models. 

Another  important  consideration  tor  rocket  exhaust  plumes  is  the  compressibility 
eflecK  I'lie  importance  of  this  effect  is  addressed  for  Case  I  (Table  7)  where  there 
is  an  infinite  velocity  ratio  between  the  jet  and  the  external  stream.  Certainly  Brown  and 
Roshko  have  pointed  out  the  importance  of  this  effect  and  the  accuracy  w  ith  which  this  effect  is 
accounted  lor  is  shown  below 

I  he  lirst  comparison  made  was  for  an  M  =2.2  air  jet  exhausting  into  still  air.  The  jet 
velocity  was  I  765  feet  second  and  the  pressure  wasambient.  I  he  density  ratio  wasasshown  in 
Table  7  tor  Case  1 

This  is  the  homogeneous  case  in  which  the  density  ratio  is  determined  by  the  Mach 
number  as  opposed  to  molecular  weight  differences  in  the  jet  and  the  external  stream.  Hence, 
compressibility  effects  are  important  for  this  case.  Since  the  kt2  turbulence  model  does  not 
account  for  compressibility  effects,  a  compressibility  correction  factor  originating  from  some 
empirical  work  tit  General  Applied  Science  Laboratory  (CiASl.)  was  utilized  to  account  for 
this  effect.  I  he  details  of  this  correction  were  presented  in  an  earlier  section. 

I'li; tire  20  compares  model  predictions  for  the  kt  2  and  k«/  turbulence  models  without 
compressibility  with  the  experimental  data.  It  should  be  noted  from  the  centerline  velocity 
profile  that  the  predicted  core  length  of  the  jet  is  too  short  compared  with  experiment.  The 
kco’  core  length  being  much  shorter  than  the  kc2  core  length.  1'his  indicates  that  the 
entrainment  of  the  ambient  stream  is  much  too  large  and  the  mixing  distance  much  too 
short  for  the  case  of  a  high  relative  velocity  between  the  two  streams. 

However,  when  the  compressibility  effects  are  accounted  for  by  utilizing  the  CiASl 
compressibility  factor  for  the  kt2  turbulence  model  and  via  the  addition  of  a  term  in  the 
modeling  equation  for  the  k«/  model,  the  results  show  a  marked  improvement.  In  fact  the 
agreement  between  prediction  and  experiment  is  excellent  lor  the  k«/ turbulence  model  with 
compressibility.  I  his  agreement  is  show  n  in  Figure  21 .  Note  that  the  core  length  prediction  is 
correct  and  the  centerline  velocity  agreement  is  excellent  out  to  60  jet  nozzle  radii.  The 
agreement  deteriorates  from  that  point  on  downstream  but  is  still  quite  reasonable.  The  kt 2 
model  with  compressibility  correction  predicts  a  core  length  that  is  too  long  and  centerline 
velocities  that  are  too  high  until  the  centerline  velocity  has  dropped  to  30  percent  of  its  original 
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Figure  20.  Mj  =  2.2  air  jet  into  stilf  air.  Comparison  of  ke2  and  ku’ 

turbulence  models  without  compressibility  centerline  velocity  profile. 


Mj  =  2.2  air  jet  into  still  air.  Comparison  of  ki  2  and  kw 
turbulence  models  with  compressibility  centerline  velocity  profile. 


v,,!uc  Hence.  both  models  agree  icasonablv  well  with  the  experimental  data  with  the  ku/ 
mi  Kiel  show  me  the  hettei  agreement  ol  the  mo  I  m  t  he  i .  it  has  been  established  that  the 
eompiessibihtx  eltccts  are  imp'inant  and  will  be  letained  loi  the  remainder  ol  tin  jet 
cotnpai  isons  in  this  section 

I  he  next  test  ot  tin  model  predictive  capabihtx  comes  about  by  comparing  the  radial 
veil  k  it  \  pi  id  lies  at  dow  list  ream  axial  stations.  1  he  1 1 1  st  station  chosen  was  near  the  end  ol  the 
ict  potential  core  as  shown  in  / tizure  2ll  at  22.9  jet  radii  downstream  ol  the  nozzle  exit.  Note 
that  the  veiocuv  protile  in  l  ii'iirc  22  shows  a  small  diameter  potential  core  at  this  distance 
Jovv nsti earn  and  a  gradual  velocitv  reduction  as  the  ladial  distance  increases.  Note  that  both 
turbulence  models  show  excellent  agreement  with  the  experimental  data  with  the  ku/  model 
showing  slightly  hettei  agreement 

I  he  next  point  chosen  loi  comparison  was  at  43.9  nozzle  radii  dow  nstream.  I  he  radial 
pml tie  at  this  axial  station  is  shown  in  Ilium'  22.  I  he  agreement  between  the  ku2  model  and 
expenment  is  near  peitect  at  this  axial  location.  The  kt2  model  prediction  shows  a  reduced 
velocitv  compared  with  experiment  indicating  that  mixing  is  occurring  slightly  too  rapidly. 

I  he  last  point  chosen  lor  comparison  was  at  61 .7  nozzle  radii  downstream.  The  radial 
profile  at  this  axial  station  is  shown  in  liiuue  24.  1  he  agreement  between  the  ku/  model  and 
experiment  is  v  cry  good  at  this  axial  location.  I  he  k<  2  model  prediction  again  shows  a  reduced 
veil).. if v  computed  with  experiment. 

It  is  tints  concluded  that  compressibility  etlects  are  important  tor  Case  1  t  Table  7i  and 
must  be  accounted  tor  to  achieve  a  reasonable  agreement  between  theory  and  experiment.  It  is 
also  concluded  that  the  ku/  turbulence  model  shows  improved  agreement  with  the 
experimental  data  compared  with  the  kr2  model 

Ihe  second  comparison  made  was  lor  a  Mj  =  0.X9  1C  jet  into  a  supersonic  M0  =  1.32 
air  jet.  The  jet  velocity  was  3520  It/scc  and  the  external  velocity  was  1295  ft/sec.  Fite 
pressure  was  ambient  and  the  initial  density  ratio  is  shown  in  Table  I  //. 

/  ign/'i  V  compares  model  predictions  tor  the  kc2  and  kw’  turbulence  models  with 
compressibility.  It  should  be  noted  that  the  agreement  between  theory  and  experiment  tor 
the  core  length  is  reasonable.  The  ku/  model  shows  a  slightly  better  agreement  with  the 
experimentally  determined  length. 

Ihe  india!  velocitv  profiles  were  First  compared  at  an  axial  distance  downstream  ol 
I  I  02  nozzle  radii.  I  his  can  be  seen  from  Tibiae  23  to  still  be  in  the  potential  core  region.  Note 
m  /  ti;iiie  2< >  that  the  predicted  radial  protile  utilizing  the  ku/  model  is  reasonable  but  certainly 
docs  not  agiee  as  closely  as  it  did  m  Case  I  above,  (discrepancies  between  experiment  and 
theory  became  apparent  when  the  velocity  diftcrcncc  ratio  has  dropped  to  approximately  60 


Mj  =  2.2  air  jet  into  still  air.  Comparison  of  ke2  and  ku'  turbulence 
models  with  compressibility  radial  profile  at  x/n  =  22.9. 


Mj  2.2  air  jet  into  still  air.  Comparison  of  kc2  and  ko>'  turbulence 
models  with  compressibility  radial  profile  at  x/n  =  43.9. 
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percent  of  its  initial  \alue.  Use  at  the  k<2  model  gives  a  noticeably  poorer  agreement  with 
experiment  over  the  entire  range  ol  comparison. 

I  he  next  radial  profile  to  he  compared  is  show  n  in  Figure  27  at  a  dow ns t ream  distance 
ol  IV  16  no//le  radii.  Note  that  the  comparisons  do  not  agree  at  the  r  =  0  point  because  the  core 
lengt  h  is  incorrectly  predicted  bv  hot  h  tui  bulcnce  models  Hot  It  models  do  equally  well  over 
t  he  entile  range  ol  para  meters  at  t  his  dow  list  ream  (.list  a  nee  w  it  h  t  he  most  serious  discrepancies 
occur!  mg  at  t  (I. 

I  he  last  radial  profile  comparison  was  made  a t  a  distance  ol  30.88  nozzle  radii 
dow  nstream  of  the  exit  plane  as  show  n  in  Figure  2i S',  lfoth  models  do  reasonably  well  over  the 
radius  range  covered  with  the  largest  diset  'pancv  occurring  a t  the  lower  values  ol  the  radius 
ratio. 


I  he  preceding  comparisons  demonstrate  that  the  turbulence  kinetic  energy  methods 
predict  tiie  velocity  profiles,  both  radially  and  axially,  reasonably  well  over  a  fairly  wide  range 
of  velocities  and  densities.  "Reasonably  well"  means  an  accuracy  within  =30  percent  tor  a 
ma ximum  error.  I  hey  also  demonstrate  that  compressibility  effects  are  important  and  must  be 
accounted  lor  in  the  modeling  in  order  to  achieve  this  accuracy  .  Otherwise,  even  larger  errors 
will  occur. 

1  hese  comparisons  also  demonstrate  that  the  k«/  turbulence  model  show  closer 
agreement  to  cypcrimcnt  than  the  02  mode)  over  the  range  ol  comparisons  made.  Hence  this 
model  can  be  considered  as  a  viable  alternative  to  the  kt2  model. 

VII.  Ri -ACTING  SHEAR  LAYER  COMPARISON 

I  hits  lar  the  turbulent  mixing  models  have  been  compared  against  flows  which  have 
ditferent  velocities  lor  e.  It  ol  the  two  streams  or  different  velocities  and  different  densities 
I  he  latter  type  flow  is  certainly  more  applicable  to  the  rocket  exhaust  plume  Hows  since  large 
density  ratios  occur  Im  these  cases  and  since  it  is  important  to  know  the  species  distribution 
across  the  mixing  layer.  I  he  next  level  ol  complexity  in  the  modeling  procedure  is  to  examine 
the  turbulent  mixing  oi  a  reacting  shcat  layer,  preferably  one  in  which  the  initial  densities  ol 
the  two  streams  are  the  same.  I  Ins  eliminates  another  variable  in  the  problem  if  the  initial 
density  i.itio  can  be  held  at  unity.  In  addition,  a  kinetically  simple  reacting  How  sys¬ 
tem  is  mandatory  for  an  accurate  test  of  the  turbulence  models  since  chemical  reaction 
tales  constantly  change. 

Reacting  shear  layer  experiments  ol  this  kind  have  recently  been  completed  at  the 
I  m v ei sit \  ol  Adelaide  in  Australia  and  have  been  re  oiled  at  the  Second  Svmposium  on 
I  inbiileiit  Shcat  Nows  bv  Wallace  and  Ifiovvn  |  S  |. 
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Figure  28.  Mj  =  0.89  H2  jet  into  M 1 .32  air.  Comparison  of  kf2  and  k<r/  turbulence 
models  with  compressibility  radial  profile  at  x/r j  30.88. 


In  these  experiments,  the  jet  stream  consisted  ol  NO  mixed  with  N.-  in  various  low 
level  concentrations  and  the  external  stream  consisted  ol  ()>  mixed  with  N;  in  various  other 
low  level  concentrations.  (lie  jet  stream  velocity  was  K2  leet/second  and  the  external 
stream  velocity  was  lo  leet'second.  The  How  channel  was  o.b4  inches  wide  and  the  height 
of  the  jet  stream  nozzle  was  ().'»N  inches  while  the  height  ol  the  external  stream  nozzle  was 
I/O  inches,  (  he  Ke  ntimher  based  on  the  boundary  layer  height  was  '  100. 

I  he  chemical  reaction  involved  in  this  experiment  is  very  simple 
NO  .  O,  -  NO  t  () 

I  his  is  a  vvcll  known  chemical  reaction  and  its  rate  is  known  with  a  high  degree  of  accuracy  as 
long  as  the  temperature  remains  below  51)0  degrees  k  I  herelore.  this  experiment  allows  the 
comparison  of  the  turbulent  mixing  model  directly  since  the  chemical  kinetics  of  this 
reaction  are  so  well  known. 

Comparisons  were  made  tor  three  two-dimensional  reacting  shear  layer  cases  shown  in 
/Vgt/rc  29  and  further  details  in  1'tihlc 

Note  in  //gi/rc  29  that  the  jet  stream  contains  the  luel  NO  and  the  external  stream 
contains  the  oxidizer  ( >>.  I  he  low  concentrations  ol  luel  and  oxidizer  mixed  with  the  carrier 


Figure  29.  Two-dimensional  reacting  shear  flow  schematic. 


TABLE  8.  INITIAL  CONDITIONS  FOR  REACTIVE  SHEAR  LAYER 
COMPARISONS 


CASE 

NO 

Uj  ue 

Pj  Po 

Tj  To 

JET  STREAM 

CONSTITUENT- 

EXTERNAL 

STREAM 

CONSTITUENT' 

NO 

N2 

03 

N2 

1 

5.0 

0.996 

1.0 

0.05 

0.95 

0.010 

0  990 

It 

5.0 

0  977 

1.0 

0.05 

0.95 

0.0379 

0.9621 

III 

5.0 

0  95 

1.0 

0.05 

0.95 

0.078 

0.922 

'  Constituents  given  as  mole  fractions 

stream.  N  is  detailed  m  /',;/>/<•  S  tor  the  eases  that  weie  compared  I  he  reason  lot  such  low 
eoneentrat  ions  ot  react  tints  is  the  degree  ot  react  i\  it\  ot  the  NO  ( ) .  reaction.  A  large  amount  ot 
heat  is  produced  trom  minute  quantities  ot  reactants. 

Comparisons  were  made  lor  the  cases  show  n  in  Tublcfi  where  theamount  of  oxidizer  in 
the  external  stream  was  constantly  increased  from  I  to  7.X  percent.  I  his  increased  the  ()  1 
ratio  and  resulted  in  increasingly  higher  shear  layer  temperatures.  C  omparisons  were  made  of 
both  the-  \  eloeity  and  temperature  profiles  lor  these  eases.  1  Jensity  profiles  were  not  measured. 

I  lie  virtual  origin  x.  was  not  determined  in  these  tests  as  it  was  lor  the  non-reacting 
shear  laxer  experiments  presented  earlier,  lienee,  for  the  theoretical  predictions.  x„  was 
taken  as  zero  indicating  that  the  shear  layer  starts  growing  exactly  at  I  he  exit  plane  of  the 
nozzle,  litis  will  produce  some  emu  In  sealing  the  width  of  the  predicted  shear  laxei 
incorrectly.  lienee,  exen  it  the  theoretical  and  experimental  icsiilts  agreed  perfectly,  the 
width  i’l  the  stieai  laxei  would  show  a  diserepa.  ..  between  the  predicted  and  measured 
results.  I  low  ex  er.  this  A  not  considered  to  tv  a  sen. ms  discrepancy  m  light  of  the  magnitude 
ot  the  disagreement  ot  tiie  temperature  profiles  as  will  he  subsequently  shown. 

Moth  the  k r  ?  turbulence  model  and  the  ko/  ttuhulenee  models  were  utilized  in  the 
predictions  ot  the  reacting  shear  layers  1  hese  models  produced  virtually  identical  results  ;ind 
therelote  the  predictions  shown  were  those  obtained  utilizing  the  km'  turbulence  model. 

M 


Results  comparing  the  reacting  shear  layer  tor  C  ase  1  given  in  Table  <3  predicted  and 
measured  are  show  n  in  Tillin’  30.  I  his  is  a  comparison  ol  the  velocity  profile  across  the  shear 
layer.  I  he  agreement  between  theory  and  experiment  is  excellent.  1  his  agreement  compares 
!a\orahl\  with  tiie  non-reacting  shear  layer  results  presented  earlier  for  the  case  where  the 
initial  density  ratio  is  I .  That  is  the  situation  here.  1  he  density  ratio  is  0.996  initially  as  shown  in 
Table  .s'.  1  herelore  it  appears  that  the  velocity  protile  is  little  attecied  by  the  reacting  flow  at 
least  at  this  level  and  similarly  good  agreement  between  experiment  and  theory  is  shown. 

l  i^itre  31  compares  the  temperature  profile  predicted  across  the  reacting  shear  layer 
will)  that  actually  measured.  In  contrast  to  the  excellent  agreement  between  the  velocity 
profile  determined  from  experiment  and  theory,  these  results  arc  quite  the  opposite.  The 
predicted  temperature  rise  profile  across  the  shear  layer  looks  nothing  like  the  measured 
profile.  1  lie  predicted  maximum  temperature  rise  is  more  than  twice  what  was  actually 
measured,  likewise  the  predicted  temperature  gradient  is  much  larger  than  that  found 
experimentally.  I  he  shear  layer  width  also  does  not  agree  but  as  was  mentioned  previously, 
this  is  scaled  In  the  virtual  origin  \„  and  therefore  cannot  agree  unless  the  origin  of  the  shear 
lay  er  lies  at  the  exit  plane  of  the  two  Hows.  In  addition,  it  is  worth  noting  that  the  location 
of  the  predicted  maximum  temperature  rise  from  the  dividing  streamline  is  displaced  toward 
the  jet  (fuel  stream)  side  relative  to  the  experimental  results. 

l  hc  results  shown  1  or  both  /Vgi/rev  30  and  31  were  for  an  axial  distance  ol  100  mm 
downstream  of  the  shear  layer  exit  plane. 

for  the  3  79  percent  ( mole)  ()>  case  given  as  Case  1 1  in  Table  S.  no  experimental  results 
were  av  ailable  w  ith  w  hich  to  make  a  comparison  of  the  velocity  profile  across  the  shear  layer. 
Hence  the  temperature  profile  is  compared  again  at  100  mm  downstream  of  the  shear  layer 
exit  plane.  I  his  comparison  is  shown  in  Tif’nre  32.  Note  that  the  essential  features  that  were 
discussed  for  the  I  percent  () *  case  shown  in  Tii’nre  31  are  applicable  here  as  well.  The 
predicted  maximum  temperature  rise  is  more  than  twice  the  measured  value.  The  location  ol 
the  maximum  is  again  shifted  toward  the  jet  stream  side  and  the  width  of  the  shear  layer  is  in 
disagreement. 

A  I  ewis  number  variation  was  made  to  determine  its  effect  on  the  width  of  the 
temperature  shcai  layer  compared  with  the  velocity  shear  layer  but  only  minor  variations 
resulted.  I  Ins  did  not  affect  the  overall  poor  comparison  between  theory  and  experiment. 

1  he  next  case  compared  is  shown  as  Case  Number  111  in  Table  <V.  The  Ov  mole  fraction 
was  increased  to  7  S  percent  lor  this  comparison.  I  he  velocity  profile  is  compared  in  T'ignre 33. 
Note  that  the  agreement,  w  hiic  not  as  good  as  tor  the  1  percent  ();  case,  is  again  excellent.  1  his 
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Velocity  profile  comparison  for  reacting  shear  layer 
Table  8,  case  number  III. 


v  lion  111  tv  expected  hovvcv fi  in  light  ol  the  results  presented  in  /  inure  29.  Again  this  is  a  nearly 
constant  imii.il  density  ease,  the  density  ratio  being  0  hg  as  given  in  Table  S. 


/ /gore  29  illustrates  the  comparison  ol  the  predieted  and  measured  temperature  profile 
aei  oss  the  reacting  shear  lay  er.  Similarly  lot  this  ease  the  predicted  maximum  tempera tute  rise 
's  nearly  twice  the  measured  value.  I  lie  predicted  maximum  is  similarly  skewed  to  the  jet 
stream  side,  ilenec  the  common  Icaiurcs  that  were  noted  in  Cases  I  &  II  (  Table  Hu  are  similarh 
ev  ident  in  (  ase  1 1 1 

It  is  instructive  to  examine  the  results  lor  an  increasing  amount  ol  oxidizer  and  noting 
the  shi  It  tow  aid  the  fuel  side  I  Ins  is  shown  in  It  wire  22  w  here  the  maximum  temperature  rise 
location  moves  toward  the  oxtdi/er  stream  monatonieally  with  increasing  oxidizer.  This  leads 
one  to  suspect  that  the  turbulence  model  utilized  in  the  predictions  is  behaving  in  a  "laminar 
manner".  Note  that  the  experimental  results  show  very  little  shif  t  in  the  maximum  temperature 
rise. 


I  his  livpothesis  was  confirmed  alter  programming  a  laminar  mixing  model  and 
introducing  it  as  an  option  into  the  IK)  A  I  analysts  1  his  was  accomplished  1 8 1  and  the  results 
are  shown  in  I  inure  2'i.  Note  that  the  laminar  mixing  case  results  in  virtually  the  same 
maximum  temperature  rise  as  lor  the  k«./  turbulence  model.  I  he  only  difference  is  the 
spreading  rate  o!  the  laminar  sheat  In  yet  compared  to  the  turbulent  shear  (aver. 

finally  predictions  were  made  with  two  other  turbulence  models;  (i)  Prundtl  mixing 
length  model  and  lii)  I )onaldson-<  iray  eddy  viscosity  mode!  1  lie  results  are  shown  in  Figures 
2  7  and  fv,  respectively .  It  is  ev  ident  1  to  in  these  ligurcs  that  neither  ol  the  simple  models  offers 
any  hop  -  ol  Ivttei  agreement  between  experiment  and  theory. 


VII!  CONC'l.l 'SIONS 

Results  ol  this  study  clearly  show  that  matching  of  a  theoretical  velocitv  profile  lot 
mean  velocities  with  experimental  results  is  not  a  good  indicator  of  the  correctness  of  a 
tni  bit  fence  model.  I  his  w  as  show  n  most  vividly  lot  the  t  eaet  mg  shear  la  vet  ex  pet  intents  vv  here 
the  velocitv  profile  match  was  excellent  and  the  temperature  profile  through  the  shear  layet 
vv  as  in  ei  i  oi  by  mote  than  a  laetoi  oltwo.  I  lie  velocity  prolile  was  the  easiest  parameter  to 
match  when  utilizing  'fie  tuibitlenee  kinetic  energy  models  given  in  this  investigation.  I  he 
constant  initial  density  cases  showed  excellent  agreement  between  experiment  and  theory  tor 
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Figure  35.  Measured  and  predicted  temperature  distribution  in 
shear  layers  between  nitric  oxide  and  ozone. 
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Figure  37.  Temperature  profile  prediction  for  reading  shear  lay n 

using  Prandtl  mixing  length  turbulence  model.  Table  8  ,i>. 
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Figure  38.  Temperature  profile  prediction  for  reacting  shear  layer 
using  Donaldson-Gray  eddy  viscosity,  turbulence  model 
Table  8,  case  number  III. 


all  cases  investigated.  I  his  was  not  true,  however,  tor  the  case  where  the  initial  density  ratios 
were  significantly  greater  than  one  (=7). 

The  density  profiles  across  the  shear  layer  were  very  poorly  predicted  for  all  cases 
examined  in  this  investigation.  I  his  was  especially  true  for  the  case  when  the  high  velocity  jet 
fluid  was  simultaneously  the  low  density  fluid  in  the  two-dimensional  shear.  This  isacasethat 
more  nearly  corresponds  to  the  rocket  exhaust  plume. 

The  k<o'  turbulence  model  gives  comparable  results  with  the  kt2  turbulence  model.  In 
all  cases  examined  the  ku/  model  performed  as  well  as  the  kt2  model  and  in  certain  instances,  it 
performed  much  better  than  the  kt2  model.  This  was  especially  true  for  the  M=2.2  air  jet 
exhausting  info  still  air. 

The  temperature  profiles  that  are  predicted  by  the  turbulence  models  are  extremely 
poor.  I  he  temperatures  are  too  high  by  approximately  a  factor  of  t  w.  This  indicates  that  there 
is  a  significant  large  structure  in  the  flows  examined.  The  turbulence  models  do  not  account  for 
this  structure  in  any  way  and  hence  are  all  deficient  in  the  basic  physics  of  the  flow.  Figure 39 
details  the  structure  of  a  typical  non-reacting  shear  layer  and  Wallace  and  Brown  |5]  have 
show  n  similar  behavior  for  the  reacting  shear  layer.  Examination  of  this  photograph  makes  it 
clear  that  the  vortex  structure  must  be  included  in  the  turbulence  models  in  order  to  obtain 
reasonable  predictions. 

.lust  how  much  large  structure  exists  in  Hows  more  typical  of  rocket  exhaust  plumes 
where  the  velocities  are  much  higher  is  not  known  a  priori.  It  is  suspected  that  the  large 
structure  will  be  less  evidenced.  If  this  is  true,  then  perhaps  the  current  turbulence  models  will 
offer  more  hope  for  making  reasonable  predictions.  However,  this  remains  to  be  seen.  Non¬ 
reactive  flow  tests  will  be  run  in  the  near  future  w  hich  w  ill  provide  the  basis  for  the  reasonable 
assessment  of  this  effect. 
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The  comparison  of  the  theory  with  the  experimental  data  of  Brown  and  Roshko[l]and 
Wallace  and  Brown  [5]  required  that  the  position  of  the  dividing  streamline  be  known.  In 
addition,  this  information  is  necessary  when  determining  the  amount  of  mass  entrained  by  the 
jet  or  by  the  shear  layer.  This  was  accomplished  in  the  prediction  program  by  utilizing  the 
subroutine  D1VSL. 

Consider  the  plane  mixing  layer  shown  in  Figure  Al.  By  taking  an  element  of  fluid 
whose  bottom  edge  is  parallel  to  the  dividing  streamline,  a  momentum  balance  along  the  fluid 
element  parallel  to  the  dividing  streamline  as  shown  in  Figure  A2  gives  the  following 

EXTERNAL  v 


JET  STREAM 

Figure  A-1.  Plane  mixing  layer. 


Figure  A-2.  Plane  mixing  layer  fluid  element  (top  half). 
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Hence 
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Hence  (A3)  becomes 


so  that 


T  x 
max 


_  *f 

pv' 2  -  pu,2)dy  +  J  Pu(uo 
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■/( 


-  u)  dy 


(A6) 


Similarly  if  we  consider  the  bottom  half  of  the  fluid  element  shown  in  Figure  A3 
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Figure  A-3.  Plane  mixing  layer  fluid  element  (bottom  half). 
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and  since 


2 

Pco  =  P  +  pv' 
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T  X  + 

max 


y*  y*  y*  _ 

/pv'  dy  =  /pu  ^u  -  )  dy  +  J pu' 

-6  -6  -6 


dy 


or 


T  X 

max 


y*  .y*  . _  _ v 

-J  pu  ^  dy  +J  p(u'2  -  v'2/ 


y* 

-6 


dy 


Equating  (A6)  and  (A9) 


y* 

f  PU  | 

(»-«.) 

dy  + 

y* 

/  p(u 

J 

-6 
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"  f  pu  ( 

\  x 7 

r 

lU  -  u  , 

V  oo  > 

)  dy  + 

J 

-6 

6 

/p( 

J 

y* 

2 

y* 

•2) 


-  v'  J dy 


dy 


and  by  assuming  V  =  v*\  (A JO)  becomes 


6  (x/  y -  ui 

/  pu(s'u)dy7pu(u-  u~  ) 

y*  (x)  2  *  X 


y*  (x) 
f  pu 
-6  (x) 


dy 


Fquation  (All)  is  utilized  to  determine  the  location  of  the  dividing  streamline. 


(A8) 


(A9) 


(A10) 


(All) 
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Figure  A4  illustrates  a  program  that  was  written  to  check  out  the  dividing  streamline 
location.  Data  from  a  realistic  2-D  shear  layer  in  the  form  of  streamwise  velocity  u and  density 
p  at  normal  locations  y  were  input  via  DATA  statements.  Utilizing  this  data  the  dividing 
streamline  was  located  utilizing  equation  (All)  which  is  coded  in  subroutine  D1VSL.  The  jet 
side  integral  11  is  given  by  the  RHS  of  equation  (All)  and  the  external  side  integral  12  is  given 
by  the  I.HS  of  equation  (All).  The  trapezoidal  rule  is  utilized  for  these  integrals. 

Detailed  output  for  the  subroutine  D1 VSL  is  given  by  the  namelist  OUT2  and  for  the 
overall  check  routine  bv  OUT1.  Note  the  location  of  the  dividing  streamline  yy  given  in  both 
namclists.This  illustrates  a  behavior  common  to  the  2-D  shear  layers  studied  -  a  bending  of  the 
dividing  streamline  toward  the  jet. 

Once  the  dividing  streamline  has  been  located,  the  entrainment  can  be  calculated  for 
both  the  jet  side  and  the  external  side.  Figure  A5  illustrates  some  additional  calculations  which 
define  two  integrals  utilized  in  the  entrainment  calculations.  The  mass  flow  integration  above 
and  below  the  dividing  streamline  is  given  by  U  and  U,  respectively.  These  integrals  are  defined 
as 


(A12) 


(AI3) 


These  integrals  are  utilized  to  calculate  the  mass  flow  changes  at  streamwise  locations 
and  differences  in  these  integrals  give  the  mass  flow  entrainment.  The  jet  side  integral  is  given 
by  U  and  the  external  side  integral  is  given  by  U. 

Detailed  output  for  these  integrals  is  given  in  namelist  OUT  3.  The  dividing  streamline 
location  is  given  in  namelist  OUTI  and  OUT2  as  before. 
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Figure  A-4.  Checkout  program  listing  for  dividing  streamline  location. 
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Figure  A-5.  Checkout  program  listing  for  dividing  streamline  location 
plus  entrainment  Integrals. 
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Figure  A-5.  (Continued). 
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Figure  A-5.  (Continued). 
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Figure  A-S.  (Concluded) . 
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For  the  reacting  shear  layer  comparison,  it  became  of  interest  to  compare  the  resulting 
temperature  rise  in  the  shear  layer  with  that  predicted  by  a  laminar  mixing  model.  Hence,  it 
was  necessary  to  add  this  capability  to  the  shear  layer  program  BOAT.  The  following 
extension  to  the  code  ( Figure  Bl)  is  necessary  to  accomplish  this.  The  details  of  the  laminar 
mixing  model  are  given  in  a  separate  report  [8]. 


117 


PRKSDXNQ  PAOS  BUMC-NOV  RLmED 


^■JfispUT  INF  FNOPTSL  (  X  I  .FX  I  »NK  ,  c  )  C)A 

C'JLATION  of  Slope  at  fjost  and  last  POINT  F0»  splice  FIT  using  lag  Clfl 
Y'wO^IALS  Cl  fl 

DIMENSION  X  I < 1 00) .FX I < 1 00)  ,S (2)  cie 


Figure  B-1.  (Continued). 
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Figure  B-1.  (Concluded). 
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SYMBOLS  (Continued) 


w 

\et  Kate  ol  Chemical  Production  ol  i  Specie 

ill 

Stream  Function  (Equation  6) 

MW 

Molecular  Weight 

k 

Turbulent  Kinetic  Energy,  '  Tu''+v'‘) 

<  1  K 

Prandtl  Number  lor  1  urhulent  Kinetic  I  nerg> 

Of 

Pin ndt  1  Number  lot  1  urhulent  Dissipation 

t  ,  r 

throughout  report) 

1  'urhulent  Dissipation  Rate 

n 

Shear  1  aver  I  hickness 

M 

Mach  Number 

a 

Sonic  Vclociu 

k 

Compressibility  Factor 

to 

Pseudo- V  orticiu 

Cp 

Heat  Capacity  at  Constant  Pressure 

u  7 

Wall  Shear  Stress  Velocity.  \  r\\ 

<5* 

Boundary  Layer  Displacement  Thickness 

IV 

Momentum  Thickness 

y 

Ratio  of  Specific  Heats 
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I  2(> 


II 


I  I  | ll;ll|'\ 


SYMBOLS  (Continued) 


kc2  Turbulence  model  constants  (Equations  10  and  11) 


kvC  ’Turbulence  model  constants  (Equations  18-20) 


Empirical  constants  in  kc2  formulation 


Eunctional  corrections  for  weak  shear  Hows  in  kc2  formulation  (Ref.  7) 


Specific  turbulent  kinetic  energy 
Defined  by  liquation  25 

Cole’s  Universal  Wave  Function  (liquation  2CM 
Wall  shear  stress 
Displacement  thickness 

r  r: 

Non-dimensional  radial  coordinate,  r - F 

re  rj 

Momentum  thickness 

Cross  stream  distance  from  dividing  streamline 
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SYMBOLS  (Concluded) 
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Of  * 
Of" 

C1 

c2 

C4 

C5 

C6 

U 


Constants  in  kto  tormulation 


Mean  axial  velocity 


V  Mean  radial  velocity 


z 

v  hIp 


max 

o 

C 

e 

j 

oo 
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Subscripts 

Maximum 
Initial  value 
Centerline 
Kxtemal  stream 
Jet  stream 
Free  stream 
Turbulent 
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